Mott transition in two-dimensional frustrated compounds 
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The phase diagrams of isotropic and anisotropic triangular lattices with local Coulomb interactions 
are evaluated within cluster dynamical mean field theory. As a result of partial geometric frustration 
in the anisotropic lattice, short range correlations are shown to give rise to reentrant behavior which 
is absent in the fully frustrated isotropic limit. The qualitative features of the phase diagrams 
including the critical temperatures are in good agreement with experimental data for the layered 
organic charge transfer salts «:-(BEDT-TTF)2Cu[N(CN)2]Cl and k-(BEDT-TTF)2Cu2(CN)3. 
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I. INTRODUCTION 

The influence of spatial quantum fluctuations on the 
nature of the Mott transition in strongly correlated sys- 
tems is currently of great interest. A class of materials in 
which these effects can be studied in detail are the lay- 
ered charge transfer salts of the k-(BEDT-TTF)2X fam- 
ily, where X denotes an inorganic monovalent anion such 
as Cu[N(CN)2]Cl or Cu2(CN)3. The electronic properties 
of these compounds have been shown to be highly sensi- 
tive functions of hydrostatic pressure fi'^'^'^'^'i'^ As a re- 
sult, the temperature versus pressure phase diagram is re- 
markably rich, exhibiting Fermi-liquid and bad-metallic 
behavior, superconductivity, as well as paramagnetic and 
antifcrromagnetic insulating phases. These observations 
suggest fascinating connections to analogous phenomena 
in various transition metal oxides.^ 

A feature of particular interest in the organic salts 
is magnetic frustration. Since the geometric struc- 
ture corresponds to an anisotropic triangular lattice, 
with incquivalent nearest neighbor hopping interac- 
tions t and t'fiSiii long-range magnetic ordering be- 
comes increasingly frustrated if the lattice is nearly 
isotropic, giving rise to an exotic spin-liquid phase 
in the absence of symmetry breakingi^^ Such a spin- 
liquid phase^iiii^ is realised in the organic insula- 
tor K-(BEDT-TTF)2Cu2(CN)3 (denoted below as k- 
CN) which corresponds to t' « 1.06i, whereas n- 
(BEDT-TTF)2Cu[N(CN)2]Cl (denoted as k-C1) with 
t' m Q.Tbt is an antifcrromagnetic (AF) insulator.-'^ 
AF order is also found in those Pd(dimt)2 salts 
for which 0.55 < t' /t < 0.85. In contrast, 
C2H5(CH3)3P[Pd(dimt)2]2 with t' = l.Obt is a valence 
bond insulator at ambient pressurCf^ Experiments on 
these kinds of two-dimensional frustrated systems have 
greatly stimulated theoretical investigations of the elec- 
tronic and magnetic properties of anisotropic triangular 
lattices is. 16.17. is. 19.20. 21. 22. 23. 24. 25. 26 

The focus of the present study is the band width con- 
trolled finite temperature phase diagram of the Hubbard 
model for isotropic and anisotropic triangular lattices. 
The key result is that small changes in the ratio t' /t 
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can give rise to fundamental changes of the phase dia- 
gram. Thus, partial and full magnetic frustration reveal 
strikingly different metal-insulator coexistence regions, in 
qualitative agreement with the experimental phase dis- 
grams for k-CI'* and k-CN^ 

The anisotropic triangular lattice has recently been 
studied also by Ohashi et alP'^ who used dynamical mean 
field theory^^ (DMFT) with a cluster extension to ac- 
count for spatial fluctuations. Although at t' « 0.8i 
reentrant behavior was found as observed for k-C1, the 
calculated Tc was much larger than the measured value. 
Moreover, only the lower boundary of the metal-insulator 
coexistence region was determined. Here, we investigate 
both the isotropic and anisotropic triangular lattices and 
use exact diagonalization'?? (ED) combined with cluster 
DMFT29 to evaluate the upper and lower phase bound- 
aries, Uci{T) and Uc2{T), of the coexistence region. As 
shown below, the shape of these boundaries, as well as 
the critical temperatures, are consistent with the exper- 
imental data for k-C1 and k-CN. 



II. THEORY AND RESULTS 

The minimal model Hamiltonian that captures the 
interplay between geometrical frustration and strong 
Coulomb interaction present in the conducting layers of 
organic salts such as k-CI and k-CN is 
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where the sum in the flrst term is limited to nearest 
neighbor sites. The hopping integrals in a unit cell con- 
sisting of three sites are ti3 = i23 = * and ti2 = t' . 
The band width is W = 9t for t' = t and W = 8M 
for t' = 0.8t. The chemical potential ^ is fixed to give 
half-filling. Within cluster DMFT the interacting lattice 
Green's function in the cluster site basis is defined as 
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FIG. 1: (Color online) Phase diagrams of Hubbard model for 
anisotropic triangular lattice {t' = 0.8t) and isotropic trian- 
gular lattice {t' = t), evaluated within ED cluster DMFT for 
t = 0.04 eV. Plotted are the first-order metal-insulator phase 
boundaries as functions of inverse local Coulomb energy U. 
In the experimental setup increasing hydrostatic pressure P 
implies increasing band width W or decreasing U. The reen- 
trant behavior found for t' = 0.8t is absent in the isotropic 
Umit t' = t. 



where k extends over the reduced Brillouin Zone and 
ujn — (2n-|-l)7rT are Matsubara frequencies. t{k) denotes 
the hopping matrix for the superlattice and E(«a;„) rep- 
resents the non-diagonal cluster self-energy matrix. This 
self-energy is calculated within ED where the environ- 
ment of the three-site cluster is replaced via a bath con- 
sisting of 6 or 9 levels, i.e., for a total cluster size Ug — 9 
or = 12. The calculations are carried out in a site 
basis and in a mixed site/molecular orbital basis. Due to 
ED finite-size effects, these treatments give results that 
differ quantitatively. Nevertheless, the qualitative fea- 
tures of the phase diagrams are consistently reproduced 
by these ED versions. Details of the cluster ED/DMFT 
formalism can be found in Refi^. 

Figure[Tfa) shows the calculated phase diagram for the 
anisotropic lattice in the region below the critical temper- 
ature. To facilitate the comparison with the experimental 
data for k-C1,* the hopping matrix elements are chosen as 
t — 0.04 eV and t' — 0.8t to reproduce the single particle 
band width, W — 0.34 eVi^i A similar value was used 
in the numerical renormalization group (NRG) DMFT 
analysis of the high-T data in Refi^. Since the data were 
plotted in a T/P phase diagram, we show the transition 
temperatures as functions of the inverse Coulomb energy. 
Increasing pressure P implies increasing electronic band 
width, so that this measurement is equivalent to keeping 



W fixed and reducing U in the calculation. The phase 
boundaries of the coexistence region are obtained by care- 
fully increasing or decreasing U from the metallic or in- 
sulating domains, respectively. Fig. 1(b) shows the phase 
diagram for the isotropic case corresponding to k-CN. 

The critical temperatures for t' = 0.8t and t' — t, 
Tc w 50 K ?s O.llt, are consistent with the measured 
values Tc w 40 K for K-Cli and « 50 K for k-CN;I 
Tc « O.li was recently obtained also for the fully un- 
frustrated square latticed On the other hand, within 
Quantum Monte Carlo (QMC) DMFT at temperatures 
T — O.lt . . A.Ot, Ohashi et ali^ found a much larger 
value, Tc w 0.3t « 140 K. The experimental data and 
the present ED/DMFT results suggest that the metal- 
insulator coexistence region is located at temperatures 
below those considered in Reff2£. 

For t' = 0.8t, the first-order phase boundaries separat- 
ing the Fermi liquid from the Mott insulator in Fig. 1(a) 
show the same kind of reentrant behavior as measured 
for K-Cl. For instance, at [/ = 1/3 eV and T « 50 K 
the system is a Mott insulator which turns into a Fermi 
liquid if T is lowered to about 20 K. Further reduction 
of T reverts the system to a Mott insulator, just as seen 
in the data. (We do not consider here the antiferromag- 
netic insulating phase which is detected at even lower 
temperature.) Ohashi et al^ found reentrant behavior 
at considerably higher temperatures. 

At present the origin of differences between the phase 
diagram for t' = 0.8t shown in Fig. 1 and the one found 
by Ohashi et al. is not clear. One reason might be 
that we consider a triangular lattice (3 sites per clus- 
ter) while in Ref^^S a square lattice with one diagonal 
was used (4 sites per cluster). However, since the ex- 
perimentally observed critical temperature is much lower 
than the range treated in Ref.^^, it would be interesting 
to apply continuous-time QMC to this problem in order 
to reach lower temperatures. 

The reentrant behavior for t' = 0.8t is in striking con- 
trast to the phase diagram obtained for the isotropic tri- 
angular lattice shown in Fig. 1(b). This limit resembles 
more closely the phase diagram derived within single-site 
DMFT.^^ The main effect of short-range fluctuations in 
the isotropic case is a significant lowering of the criti- 
cal Coulomb energy. Here, Uc2 ~ 1/2.63 eV « 9.5t, 
whereas Uc2 ~ 12t...l5i in local DMFT for the tri- 
angular lattice. Comparing Figs. 1(a) and (b), it 
is evident that anisotropy causes a further lowering of 
the critical Coulomb energies. This trend is consistent 
with Uc ~ 6t for the fully unfrustrated square lattice^i^^ 
which is topologically equivalent to the triangular lattice 
in the limit t' = 0. 

It is interesting also to analyze the width of the metal- 
insulator coexistence region obtained by increasing vs. 
decreasing pressure. For k-C1, it is observed at P « 
200 .. . 400 bar, which according to the high-T NRG anal- 
ysis corresponds to band width changes of about 2% A 
The calculated coexistence region shown in Fig. 1(a) is 
only slightly wider than this experimental range. 
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FIG. 2: (Color online) Nearest-neighbor spin correlations in 
isotropic and anisotropic triangular latticees for T — O.Obt, 
t' = ti2, t = ti3 = t23. Strong enhancement of spin correla- 
tions occurs for moderate deviations from the isotropic limit. 



The qualitative change in the phase diagram caused by 
reduced geometrical frustration can be understood by an- 
alyzing the magnetic properties of the frustrated lattice. 
In the [/ — > cx) limit, the Hubbard model can be mapped 
onto the anisotropic Heisenberg model with J' /J = 0.64 
and J'/ J = 1 for t' = 0.8t and t' = t, respectively. At 
T ~ 0, t' ^ t yields long range antiferromagnetic (AF) 
order of the 120" type, whereas t' = 0.8t gives rise to row- 
wise AF Neel order. '^^ However, in the Heisenberg model 
the temperature scale for AF order in the isotropic trian- 
gular lattice is expected to be strongly suppressed relative 
to the square lattice. For the Hubbard model, the cluster 
DMFT provides information on how the magnetic corre- 
lations < SizSjz > vary across the Mott transition in 
the isotropic case compared with t' = 0.8t. The results 
shown in Fig. [2] demonstrate that spin correlations are 
strongly enhanced as the geometrical frustration is sup- 
pressed. The isotropic lattice displays weak AF coupling 
for any U. This is in contrast to t' = 0.8i, for which the 
weaker hopping amplitude displays ferromagnetic corre- 
lations whereas spins with the larger hopping amplitude 
are antiferromagnetically coupled, indicating a row-wise 
AF Neel arrangement of spins. Thus, t' = 0.8t induces 
a much stronger tendency towards magnetic order than 
t' = t, which explains why the reentrant behavior occurs 
for t' = O.St, but not for t' = t (see Fig. [J). At low T, the 
electron entropy is suppressed for t' — 0.8t as compared 
to t' = t. As T is increased for t' = 0.8i, the system 
lowers its free energy by transforming to a metal since 
the entropy of the metal exceeds that of the ordered in- 
sulator. At even higher temperatures the system gains 
entropy of log(2) by transforming back into a paramag- 
netic insulator. This result is analogous to the one found 
for the unfrustrated square lattice. "^^ In the isotropic lat- 
tice magnetic ordering is suppressed and the reentrant 
behavior disappears. 

To illustrate the first-order nature of the metal- 
insulator transition we show in the upper panel of Fig. [3] 
the spectral weights of the cluster sites at Ep = as 
functions of U . The lower panel shows the average double 
occupancy docc = Si(^iT'^u)/3- Both quantities exhibit 
hysteresis for increasing and decreasing U , indicating co- 
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FIG. 3: (Color online) Hysteresis behavior of spectral weights 
AiiQ) of cluster sites at Ef ~ and average double occupancy 
docc as functions of Coulomb energy for anisotropic triangular 
lattice {t = 1, t' = 0.8t, T = 0.05t). Red (blue) curves; 
increasing (decreasing) U. 



existence of metallic and insulating solutions. 

Finally, Fig. [4] shows the spectral densities at Coulomb 





FIG. 4: (Color online) Average spectral distributions of 
cluster sites below and above the Mott transition for the 
anisotropic triangular lattice {t — 1, t' = 0.8t, T — 0.02t) 
for Coulomb energies U = 7.5t and U = 8.2t. The bare den- 
sity of states is shown by the dashed blue curve. 
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energies below and above the Mott transition for T = 
0.02i and t' = QM. Pfotted is the average over the three 
inequivalent sites within the unit cell. Since we are here 
concerned with the metal-insulator transition we give the 
ED cluster spectra which can be evaluated without re- 
quiring extrapolation from Matsubara to real fr'equen- 
cies. In the metallic phase the spectra reveal large quasi- 
particle weight at low frequencies as well as upper and 
lower Hubbard bands at high frequencies. The insulating 
phase exhibits a Mott gap, as well as pronounced spectral 
weight in the region of the Hubbard bands. Qualitatively 
similar features are also seen for the unfrustrated square 
latticci^i^i 

III. CONCLUSION 

In conclusion, the phase diagrams of the Hubbard 
model for the isotropic and anisotropic triangular lat- 



tices have been determined within cluster DMFT and ex- 
act diagonalization. For moderate frustration, t' = 0.8t, 
reentrant behavior is found and the phase boundaries of 
the metal-insulator coexistence region are in qualitative 
agreement with the T / P phase diagram observed experi- 
mentally for the anisotropic organic salt k-C1. The reen- 
trant behavior disappears in the fully frustrated limit, 
t' = t, in agreement with measurements on the nearly 
isotropic compound k-CN. The phase diagram then bears 
overall resemblance to the one obtained within local 
DMFT, i.e., in the absence of inter-site correlations. The 
critical temperatures, Tc « 50 K for the isotropic and 
anisotropic lattices, are consistent with the data for n- 
CN and k-C\. These results should also be relevant for 
the phase diagram of [Pd(dimt)2]2 salts exhibiting small 
deviations from the isotropic lattice.'^* 
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